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Abstract 


The altitude dependencies of the moduli of the electric field E in 
the VLF and LF frequency bands (Fb < F < /#) and in the altitude 
range of the ionosphere Z = (400 to 2500) km up to Z = 6000 km 
(the bottom of the magnetosphere) were calculated by the linear the- 
ory. The amplitudes of the field have large maxima in four regions: 
the Axis field |£o| close to the direction of the Earth’s magnetic field 
line Bo, (3 ~ 0 degrees, the fields |£s*|, and in the 

Storey, Reversed Storey and Resonance cones, /3 ~ (0 — > 20) degrees . 
Their maxima are very pronounced close to the low hybrid frequency 
Fl, Especially large one is the enhancement of the Axis field jjEbl at 
Z > (1000 — 2500) km. It grows up for about 10 6 times in this region 
and it is concentrated in very small angles A/? < 1 degree . The spatial 
distribution of \Eq\ becomes similar to a laser beam. 

The nonlinear heating of a magnetoplasma under the action of an 
electric field Ee lu>i is recently expanded by the microscopic theory by 
the author. The velocities, collision frequencies and temperatures of 
all the constituents of a magnetoplasma - electrons, ions and neutral 
particles - are taken into account. Formulae and numerical results, 
particularly obtained by a computer, are presented for the ionosphere 
in the frequency band F = (1 to 10 4 ) kHz and altitude range Z ~ 
(100- 1000) km. The temperatures are growing up very quickly with 
altitude and reach at Z ~ (150 - 200) km the ionization potential 
even when the electric field E is about some units of mvjm . At 
these altitudes T e > 10 2 • T h q. Some results of calculations by the self 
consistent solution of the basic system of equations are also discussed. 
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I. Introduction 


The proposal P2159-8-89 for the grant NAG5-1340 mainly envisaged the- 
oretical study: 

of the structure of the electromagnetic field in the ionosphere at altitudes 
(500 - 2500) km in the frequency band F = (1 - 16) kHz generated by the 
PVP Transmitter of the Sub- satellite of the Active System; 

and of the heating (temperatures) of the surrounding magnetoplasma 
under the action of the radiation of the powerful VLF generator, F ~ 10 kHz 
of the Active Satellite. 

Besides it was supposed that the study of this two problems will be done 
for realistic models of the ionosphere, based on measurements performed by 
the Active System. So, the main task of this proposal was theoretical treat- 
ment of experimental data. The proposal was for 2 years for the period 
of October 1, 1989 through September 1991 and was officially approved by 
NASA in December 1989. 

However, because of the failure and delayed launch of the Active Satel- 
lite and of the Sub- satellite by the USSR and Czechoslovakia, the term of 
this proposal was shortened to the period of one year. The author had to 
change the programme of the theoretical calculations. In general, the study 
of this problem in more detail remains very important and interesting both 
for this field and because these results can and will be used in the future by 
more successful launches of a system of two satellites with an improved pro- 
gramme, similar to the Active mission, what should be done by NASA indeed. 

Following to the main task of the proposal, the study of the structure 
of the electric field was extended for a wider band of frequencies F ~ 
(1 to 30 — 40) kHz , and a wider range of altitudes up to Z ~ 6000 km. The 
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study of the structure of the electric field at higher altitude is very important 
because the trajectory of propagation of ELF, VLF and LF electromagnetic 
waves in the magnetosphere are passing close to the Earth’s magnetic field 
lines. The author believes to extend these calculations up to the apogee of 
the magnetic field lines (tens of 1000 km) in the next study. Besides, we re- 
stricted ourselves to learn in detail only the, so called, ’’Axis field” | £^o| •> be. 
close to the direction of the magnetic field Bo- This is the most important 
and new peculiarity of the problem learned by Alpert, Budden and Moiseyev 
(see [2] to [5]). Namely the theoretical basis of these papers was used in this 
report. These calculations were done by the computer of the Goddard Flight 
Center in collaboration with James Green and Lara Aist-Sagara. 

The recently extended by the author microscopic theory was used for the 
calculations of the heating under the action of the electric field of the Magne- 
toplasma - of the temperatures of the ionosphere. By theoretical investigation 
of the parametric and nonlinear phenomena, which could be produced in the 
Tethered Magnetospheric Cloud (see [6]), the author had specially widened 
the theory, keeping in mind the interest to learn these phenomena in a broader 
altitude range than the altitude of the Tether System Z = 300 km (see also 
[7], [8]). The main part of the presented here results of calculations were 
done by the computer and program composed by Robert Estes. They are 
given in the following sections. 
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II. Altitude and frequency dependencies of the electric 
field \E\ of e.m. waves in the ionosphere Z = (400 to 6000) km. 


The altitude dependences of the moduli of electrical fields \E\ were cal- 
culated for the model of the inosphere given in Table I. This model is nearly 
conformed to the conditions of the day-time middle altitude ionosphere. The 
notations used in the table and below in the paper are: 

Z is the altitude above the Earth’s surface, N n and N e ,i are the densities 
of the neutral particles, electrons and ions; 

T°I\ and B, 7 are the temperature in Kelvin degrees and the magnetic 
field in gammas 7 (nT) along the Earth’s magnetic field line; 

v en and v e i are the collision frequencies between the electrons and the 
neutral particles and the ions; 

fo and F 0 are the Lengmuir plasma frequencies of the electrons and ions; 
f B , F b , Fl are the electron’s and ion’s gyrofrequencies and the low hybrid 


frequency of the ions 


/ 


\ 


F l = 




IbFb 

1 + f / 


1/2 


(1) 


Up to the altitude about (800 - 1000) km the ionosphere consists of 
a noticeable number of Oxygen and Helium ions (0+ and He + ) including 
protons Hi . Therefore, in the lower part of the discussed region of the 
ionosphere, for calculations of the electron and ion gyrofrequencies Jb and 
Fb , the effective mass of the neutral particles is used: 


M e fj = Mh 1 



Mhi 

Mi, , 


-1 


s = 1,2,3. 


( 2 ) 


In (2), Ni 3 and M ia are, respectively, the ion densities and masses of Hf, 
He + and Of- 
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By the numerical calculations, presented in Section II.3, the altitude • v 
rameters of the magneroplasma given in Table I are used. In Section i. J, 
where a general view of all the cones and branches of the electric field is 
given, the conditional model of the ionosphere is the following: 


So 

^ei 


3.39 • 10 6 Hz, f H = 1-194 • 10 6 Hz f L = 2.49 • 10 4 Hz, 

10 -2 sec -1 , ^ = — — = 5.44 • 10 -4 (3) 

M Hl MeJJ 


II. 1 Statement of the problem. Equations, formulae. 

The linear theory of radiation of an electric dipole in a magnetoplasma, 
used in this study, was developed by Alpert, Moiseyev (1980 a,b, 1983), and 
Alpert, Budden, et al. (1983) (see [2] to [5]). The theoretical results of these 
papers, namely: the main formulae, the computer program and the general 
physical understanding of this problem, obtained in these studies, are used 
here for producing a picture of the expected structure of the electromagnetic 
field radiated in the neighborhood of a mother satellite, moving in the iono- 
sphere at altitudes of the Active system. Besides it is supposed that this 
field is recorded on the far 2one of this source by a child sub-satellite, moving 
around the mother satellite at distances from it r ~ 10s to 100 km or a 
little more. Then the developed in the cited above papers theory for a ho- 
mogeneous medium may be used for calculations of the electric field around 
the mother body, at the discussed here altitudes of the ionosphere and in 
the magnetosphere . This is a sufficiently good approximation because at 
the interesting for us altitudes, on such distances r from the mother body, 
the magnetosphere may be considered as a homogeneous medium. It is also 
supposed that the local parameters of the magnetoplasma (the electron and 
ion densities, the collision frequencies, the magnetic field, etc.) in the neigh- 
borhood of these bodies were estimated by this experiment simultaneously 
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with the recording of the amplitudes of the electric field. 


An homogeneous cold non-magnetic magnetoplasma, which is character- 
ized by a superimposed magnetic field, and by an electrical tensor £ 0 and a 
refractive index n is considered. The magnetic field is parallel to the z axis 
in the cartesian ( x,y,z ) and in the cylindrical (p, <p, z) coordinate systems. 
The refractive index may be thought by us as a vector n with components 


n x = nsin0cos<p, n y = nsinO • simp, n z = n cos 0, 
K = (n 2 x + nl) = n 2 sin 2 Q, 


( 4 ) 


where 0 is the angle of the wave normal to the vector k with the z axis, 
z || Bo- The electromagnetic waves E,H ~ e ,wt , generated in the magneto- 
plasma, are produced by an electric dipole of moment Ie ,u,t , uj = 2nF is the 
angular frequency of the waves. The electric dipole is parallel to the z axis, 
i.e. 1 1| Bo- The source dipole is placed at x = y = z = 0, i.e. on the mother 
satellite. The receiving point is at a distance r = y/ar* + z 2 from the source 
in a direction ft to the magnetic field Bq and 


x = r sin ft, y = 0, z = rc os ft 
since the plasma has rotational symmetry around the z axis. 


( 5 ) 


The general solution of the system of the Maxwell equations of this prob- 
lem is described by a sum of complicated integrals (see [2], [3]). The inte- 
grands of these integrals are rapidly oscillating functions. They are described 
by the Bessel functions J 0 *p) and J\ and the derivatives 

where 

nx = rcsin0, n.|| = ?icos0, n 2 = n 2 ± + njj. 

These integrals may be studied only by numerical methods. The method 
of the steepest descents, i.e. the stationary phase method was used for this 
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purpose. Certainly, the accuracy of this asymptotic method is sufficiently 
high only in the far zone from the source, namely, when 


— \fp 2 + r ^ 1- (6) 

c v c 

The main contribution to the field is made by saddle points. Two cases 
should be considered by analyzing the integrals. Namely these cases charac- 
terize the basic physical properties of the field. 


The first case is when the observation point - the receiver - is at very 
small horizontal distances from the source, namely, when x ~ 0, i.e. the 
angle 0 of the ray direction is very small. In this region the field is growing 
up and becomes very strong close to the direction of the magnetic field. It 
is called the Axis field. The saddle points are estimated in this case by the 
equation 

= 0 . (7) 


dn || 
dn i 


The Axis field enhancement appears when u > u?l, ul, is the low hybrid 
frequency. This is one of the most interesting peculiarities of this problem 
(see [1] and below). For many decades, an erroneous conclusion was in the 
literature that the field disappears close to the direction of the magnetic field 
B 0 (see, for example, Arbel and Felson [5]). 


The second important case is when n x p 1 and the observation point 
is sufficiently far from the axis, from the direction of the magnetic field Bo- 
In this case, the asymptotic of the Bessel function 

Jo, Ji ~ (2x n x p) 1/2 exp -i (ii x p - (8) 

may be used. The saddle points are estimated in this region by the equation 


dn 

cos 0~r~^~ + sin 0 = 0 
dn j_ 


( 9 ) 
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The field is enhanced in the regions where two saddle points of the inte- 
grands are close to each other - they coalesce. Two regions of enhancement 
of the field exist in this case in different angle intervals /?. They appear in 
the, so called, Storey and Reversed Storey cones. The field becomes espe- 
cially strong in the Reversed Storey cone at frequencies w < ui (see [2] to [5]). 

Certainly, one more region of the enhancement of the field is the resonance 
cone. It occurs when the coefficient of refraction n — » 0. At the resonance 

0 — /? = db — , tan 2 /3 = — (10) 

“ £zz 

where £ xx and £ zz are the elements of the tensor £o- In our case, the cartesian 
coordinates 

£ j-j ^ xy 0 

£() ^yx £ yy 0 ? &xx ^ yy (11) 

0 0 £ zz _ 

The computer program used in the cited above paper by Alpert, Budden, 
etc. [3] for numerical calculations of the moduli \E\ of the electric field in 
the four regions of their enhancement, was developed by K. Budden, the for- 
mulae for the components of the electric field given in this paper were used. 
They remain complicated even by using of the first order descents evaluation 
of the integrals which gives the contribution to one saddle point of the Axis 
field [E’ol and to two saddle points of the Storey |i?st| and Reversed Storey 
cones. 

The field Eq is expressed by an algebraic combination of the both Bessel 
functions and Ji( n ±iP) (see (5)), and also by a combination of the 

components of the refraction index n and of the elements of the tensor e of 
the plasma. The fields I and l-EflevStl are expressed by a combination of 
the Airy integral function Ai(p) and its derivative Aj(p), and also by n and 
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e 0 . The field in the resonance cone |£^fl es | is much simpler and is expressed 
by an algebraic combination of the elements of the tensor, and the refraction 
index of the plasma. Certainly, all these formulae depend on the angular 
frequencies u and the angle of the ray direction /?. 

For orientation the formulae used by the programming are given here 
only schematically. They give a general presentation about the distance, 
frequency and angle dependence of the interesting for us moduli of the field 

\E\ = (\E x \ 2 + \E y \ 2 + \E z \ 2 ) 112 (12) 

These asymptotic formulae are of the following shape: 


The Axis field 



• E Ax [J 0 , Ji, n, c 0 ], ( 13 ) 

where all the values of F Ax \- • ■] depend on u, /3 and on the characteristics 
frequencies of the magnetoplasma (see [2]). 


The field in the two Storey cones 

u; 3 \ y/n ■ exp(—i {. . .}) 


Est * / 


' ^ ^ (“r) ^ • (cos fl) 1 / 3 ■ (sin/?) 1 / 2 

Fst [a,(p), A'(p), n, co] , 


(14) 


where all the values of Fst[- • •] depend on u>, ft, etc., Ai(p) and A’(p) are the 
Airy function and its derivative, 



and {...} is determined by some parameters of the magnetoplasma, by u 
and 0 (see [3]). 


The field in the Resonance cone 


I-Eflesl — M ^ 


n p n z exp — i-r(n p s\n0 + n z cos /?)j 

' ~W) 


• FR ea [n(w,/?), e o (u,0), 0] . (16) 

Let us note here the sufficiently different dependence of |.£'o|, l-^stl, l-EflevStl 
and on distance from the source. 


II. 2 Frequency dependence of \E\ in all regions-cones of their 
enhancement. General picture. 

The distribution of the cones of the electrical field amplitude \E\ enhance- 
ment, formed in the magnetoplasma around the direction of the magnetic 
field Bo, i.e. around the Axis field |Eo| in all the frequency bands is shown 
on Fig.l. 

The generatix of the Storey cones containes the angles 

0M = (0St,M t° 0St,o), (17) 

where 0st,M is the maximal value of 0m of the maxima of the angle depen- 
dence /?(0), i.e of the direction of the group velocity f/ = — of the wave. 
Let us remind that 0 is the angle between the wave vector k and the mag- 
netic field Bo- By both of these angles the 0(0) dependences have points of 
inflection. The generatix of the Reversed Storey cone contains the angles of 
the minima of the dependence 0(0), they are 

0m = (0St,m to 0). (18) 
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The frequency bands of these branches of waves are Fm = (Fst,M to fso) 
for the Storey cone and F m = ( Fst,M to Fi) for the Reversed Storey cone. At 
the frequency FstMi where Fb < Fst,M Fl, the Storey cones are combined 
(see below Fig. 4 and [3], [4]). The frequency Fso lies in the band of the 
electron whistler mode, fso < /s- For the following in this section results 
of calculation of the fields \E S t\ and |F nevS t|, these characteristic frequencies 
and angle are equal to: 

FstM = 2.6 ■ 10 -3 /s, 0st,M ~ 21°, 

F S( ,o = 0.4 S/b, F b = 5.4 • 10- 4 /b (19) 

and Fb is the ion’s gyrofrequency. The most important characteristic of the 
Storey cones fields is the large enhancement of IF/ieaStl i n th e Reversed Storey 
cone by approaching the low hybrid frequency, namely, when 

Fst,m —> Fi, f3 m — + 0. (20) 

The resonance cones are formed in the angle range /3 Bea = (0 to 90)° in the 
frequency band Fi to f B - 

The frequencies used in this study lie in the very low frequency band 
VLF, F B < F < Fi and in the low frequency band LF, F B < f < f B . By 
calculating of \E\ in the extra low frequency ELF, 0 < F < F B , the pres- 
ence of ions of different kind should be taken into account. In general, the 
behavior of \E\ of ELF waves is the same, however, it becomes more compli- 
cated because of the presence of several ion gj'rofrequencies, ion-ion hybrid 
frequencies and crossover frequencies (see [4], [5]). The angle and frequency 
dependencies of the moduli of the electric field |F| of VLF and LF waves 
presented here for the model ionosphere, given above by (3), are shown on 
Figs. (2 to 4). Only two regions of enhancement of the field exist in the VLF 
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band. They are the branches l-EfleuS*! and \Est\- Depending on /?, they are 
oscillating when F = const and are given on Fig. 2 for F = 10 -2 /b < Fl- 
Three branches |£o|> |-Est| and \Efu:,\ exist in the LF band. The field |£o| is 
also an oscillating function on /?. Three fields are given on the same figure 
for Fl < F = 10 ~ 2 Jh- The angle dependence of |£/?e S | for different fre- 
quencies is shown on Fig. 3 in more detail. It is seen that the values of the 
maxima \ER e3 \ m ax are changing sufficiently slow with frequency. However, by 
approaching the low hybrid frequency, the value of \E Res \ max , similar to the 
\E Rev st\max and \E 0 \ max , is growing up very quickly: \E 0 \ max becomes about 
10 3 - 10 5 times larger by small changes of the frequency close to F L . The 
shape of the frequency amplitudes of the electric field dependencies in this 
region is of resonance type. This is seen in more detail in the next section, 
where results of calculation of the Axis field Eq are given. 

Frequency dependencies of the maximal values of the moduli of the elec- 
tric field \E\ max in all the regions of their enhancement are shown on Fig.4. 
They are calculated for the model of the ionosphere given by (3). 


II. 3 Electric field |£ 0 | nearby the Axis, direction z || B. 

The behavior of the electric field around the source of radiation at dif- 
ferent altitudes of the ionosphere is presented in this section by detailed 
calculation of the Axis field IEqI. These results are presented here by the 
Figs. (5 to 10) and Table II. By calculation of \E\ in both this and previous 
section, the factor (/^) is omitted (see formulae (5) to (9)). This factor is 
used only on one of the dependencies shown on Fig.10. 

The oscillating character of |£ 0 |, depending on the horizontal distance x 


15 



4 


5 6 7 8 9 K) 11 12 
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Fig 2. Angle dependencies of the moduli of the electric field |i?o|, \EllevSt\ 
and | Efes\ along the axis, in the Storey and Resonance cones. 
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Fig 3. Angle dependencies of the electtic field moduli \Ene t \ in the Reso- 
nance cones at different frequencies far from fi. 












from the Axis z, i.e. from the direction of the magnetic field Bo at a 
distance r, ( xjr = tan 0) from the source is shown on Fig.5. Two important 
peculiarities of the field |J?o| are seen at once by looking at these plots: 

1. The amplitude is growing up very quickly especially in the altitude 
region Z = (500 to 800) km. 


2. The maximal values of \Eo\ max are very close to the direction of Bo, 
x < 10 -1 km, 0 ~ * <C 0.1°. 

The frequency (F/Jb) dependencies of Max\Eo\ max , i.e. the maximal 
values of the maxima of \Eo\ max (see Fig.5) at different altitudes are shown 
on Figs. 6 and 7. The asymmetric resonance-like, but of spiking shape of 
\Eo\ max is especially very pronounced at Z > 1000 km. By Z = 2500 km 
and Z = 6000 km, the values of |£ 0 | are changing (10 3 - 10 4 ) times in very 
small intervals. In general, the amplitude at these altitudes is diminished for 
about (10 s — 10 6 ) times up to the more or less stabilized values of |.Fo I max- 
The angle dependencies 0(F) are illustrated on Figs. 8 and 9, where the 
l^olmax dependencies are also shown. The spiking shape of 0(F) is similar 
to the shown of the frequency dependencies of |i?o|max- The minimal values 
of 0(F) are very small, 0 <C 10 -1 degrees. 


The result of calculations, presented on Figs. (6 to 9), are summarized 
on Fig. 10 and Table II, where the altitude dependencies of A/a:r|-Elo|max an< ^ 
of its normalized value to Z = 500 km 

(Max\Eo\ max ■ F max )z 


& | Fq | m ax — 


( 21 ) 


(M a.T|£’o| max Fyyi q X ) Z— SOO hm 
are given. The frequencies F max of the corresponding values of Max\E 0 \ ma x 
are also given in the table. 
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Table II. 


Z ,km 

Max|Eo |max 


£|Eo |max 

400 

6.59 • 10" 2 

38.430 

0.08 

500 

2.920 

25.081 

1 

800 

8.53 • 10 2 

22.601 

2.21 - 10 2 

1000 

1.377- 10 5 

20.713 

2.66 • 10“ 

1500 

2.171 • 10 6 

16.588 

2.15 • 10 5 

2000 

1.583- 10 7 

13.019 

-4 

GO 

♦— * 
o 

cn 

2500 

7.160- 10 7 

10.334 

1.71 • 10 6 

6000 

1.671 • 10 8 

4.278 

2.83 • 10 5 
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|E 0 | Z = 500 km, f = 32.27 kHz 




0.00 0.01 0.02 0.03 

X In km 



0.00 0.01 0.02 0.03 0 04 


X In km 



0 00 o 02 0.04 0.06 0.08 0 '0 

X In km 


Fig 5. Dependencies of the Axis yield moduli |£o| on distance x from the 
magnetic field line Bo at different altitudes Z of the ionosphere at fixed 
frequencies. 
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)Eo|max 


Z = 500 km, r = 100 km 



(Eojmax Z ~ 800 km, r = 100 km 



F/f B 


Fig 6. Frequency dependencies of the maximal values of 
\E 0 \ max (see Fig.5) at Z = (800, 1000) km. 


99 




10 .^c 1 1 1 1 

0.0221 0.0222 0.0223 


F/f„ 

Fig 7. The same as on Fig. 6 at Z — (1000, 1500, 6000) km. 
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|Eo|max 


Z = 2500 km, r = 100 km 




Fig 9. The same as on Fig. 8 for Z = 2500 km. 


25 


Fig. 10 



Z in km 


Fig 10. Altitude dependencies of the maximal values Max\E 0 \ max and of 
their normalized values 8{M aa^Eo^^) to Z = 500 km (see (21)). 
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III. Nonlinear heating of a magnetoplasma 
by an electric field Ee ,ul< . 

It was noted above that the nonlinear theory of a magnetoplasma heating 
under the action of an electric field Ee 1 " 1 was recently extended by the au- 
thor (see Alpert, 1990 [6], 1991 [8]). In connection with the Active mission, 
theoretical calculations were performed on the basis of this theory for a re- 
alistic altitude model of the ionosphere and in which includes frequencies of 
the VLF and LF bands F = (1 - 16) kHz of this mission. The appropriate 
results of the numerical calculations are presented in this section. The ana- 
lytical formulae used by some of these calculations are given below. However, 
a sufficiently full study of this problem is possible only by a computer, be- 
cause of the big completeness of the solution of two interconnected systems 
of linear equations. For this purpose, a special computer programme was 
developed by Robert Estes, and part of numerical calculations were done in 
collaboration with him. 


III.l Microscopic theory, formulae. Limits of its applicability. 

In the microscopic (hydrodynamic) approximation, to calculate the tem- 
peratures of a magnetoplasma, containing three kind of particles: electrons 
< e >, ions < i > and one kind of neutral particles < n >, the following 
two interconnected linear systems of differential equations of the first order 
should be solved: 

Hi = _ — _ — (V e X B) - ^ e i(V e - V.) - J/ en (V e - V n ), 
dt m me 

= ^ + ^(v i xB) + ^(v.-v i )-MV-v„), 
dt M Me M 

MN n H^ = mN e v en {V e - V n ) + MNii/i n (Vi - V n ), (22) 

dt 
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and 


dT e 2 

-jit = -^eVe • E - S ei v ei (T e - Ti) - 8 en v en {T e - T no ) 

= |eVi • E + 8 ti u ei (T e - Ti) - 8 in v in (Ti - T no ) 

dT n 

= N e d en v en (T e - T no ) + NiV in 8 in {Ti - T no ) (23) 


The following notation are used in (1) and (2): 

V e , V, and V n , T e , T, and T n , being as functions of E and u> are the 
velocities and temperatures of all the kind of particles; 

N e , Ni and N n , and m, M t = M n — M are the densities of these particles; 
den-, dei, d{ n are, respectively, the energy lost by the electrons due to their 
chaotic collisions with the neutral particles and ions and energy lost by the 
ions due to their collisions with the neutral particles (the values d en ^ d t i ~ 
2 • 10 -3 and di n ~ 1 were used). 

The collision frequencies between all the kinds of particles are: 


v e ,(E,u) 


Ven{E,u) 

Vi n {E,U>) 


V e i,0 


^en, 0 




Te 0 

T e (E), 

X(E) 

i T e fl 

Tj{ E) ' 

T ifi t 


3/2 


1/2 

I 

1/2 


In 

[0-37|fL] 

In 

A 07 T’e.O 

[ 0 


(24) 


Besides, before the action of the field, at the moment t = 0, the magneto- 
plasma is isothermal T e o = 7)o = T n(Jl and quasi neutral N e0 = iV,' 0 . 

The analytical formulae of V of the solution of the system of equations 
(21) - (23) are transparent only when the dependence of the collision frequen- 
cies on the electric field is not taken into account - i.e. by the approximation 
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v = const. These formulae are given in [6], [8] when = 0, i.e. when the 
third equation of (2) is not taking into account, T n = T„o. They become 
much simpler when the angle 0 between the electric field E and magnetic 
field B 0 is zero and E || Bo- The following coordinate system is used here: 

B x = By = 0 Bo || z, 

E x = 0, E y = Esin 0, E Z = E cos© (25) 

Nevertheless, these analytical formulae are still sufficiently complicated when 
the angle 0^0. Therefore, numerical calculations with those formulae were 
also done by the computer. By learning of these formulae it was discovered 
that the used in the literature formulae (see, for examples, Gurevich, 1978 
[10]), when Bo ^ 0 can be used only when 

u 2 > u>l + v e v in , (26) 

where u denotes the angular frequencies. 


By the computer calculations, it was found that the influence of the mag- 
netic field becomes sufficiently effective only by 0 > 30° (see below). There- 
fore, the formulae for an isotropic plasma (Bo = 0, 0 = 0) may be often used 
for evaluation of the velocities and temperatures T(E). Namely 
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where i/ e = v t { + ^in and n — Again, it was found that the formula 


F e = R e ( V e ) = 
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v r 


m u> 2 + v 2 

used in the literature (see [5]) is effective only when 

W 2 > 0'e< + Vin)Vtn 


(30) 


(31) 


The formulae of the temperatures of all the particles consist of two parts. 
One is periodically changing in time at the double frequency 2u of the elec- 
tric field Ee lw( . The second part is stationary and is much larger than the 
depending on time part. The stationary temperatures T(E ) are of the most 
interest in this study. In an isotropic plasma (0 = 0), the stationary tem- 
peratures T e and 7, are estimated in the approximation 7 = const , = 0, 
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where the squared characteristic field of this problem 
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By comparison of (32) and (33) we can see that the temperature Ti(Eo,u > ) of 
the ions is smaller than the temperature T e (E 0 , w) of the electrons. However, 
really in the ionosphere this is remarkable only at the altitudes Z < (100 — ► 
200) km. At Z > 200 km, Ti(Eo,w) ~ T e (Eo,uj). 

The formulae (32) to (36) were used for some calculations of the tem- 
peratures. However, an important part of these calculations was done by 
computer programme. This helped us to learn the limits of the suitability - 
applicability of the approximations: 

v = const , = 0, T n = T n0 , 

at 

v — i/(E,w), -77 ^ 0, T n = T n (E,u;). (37) 

at 

These results are given in the next section. However, first of all, let us single 
out here what, in general, are the limits of applicability of the microscopic 
hydrodynamic theory of the heating of a magnetoplasma under the action of 
an periodically alternative electric field, used widely in the literature. 

It is obvious, that in a collisional medium, without any other losses, the 
system of the equations (2) of the derivatives ^ and has not a 
stationary solution. Indeed, its determinant A3 = 0. Even the temperature 
of the neutral particles is growing up continuously in time under the action 
of the electric field and becomes a source of heating of electrons and ions 
(see below Fig. 17). By shortening of this system, namely, assuming that 
^ = 0, T n = T n0 , the solution of the shortened system of two equations 
becomes stationary because of additional source of losses of this process be- 
comes the saved in the equations collisions v en and i/,- n between the electrons 
and ions with the neutral particles. On Fig. 11, the establishment in time 
of the temperature T e (E 0 ,u>), calculated by the selfconsistent solution of the 
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full system of equations (22), and of the shortened system of equations (23) 
are shown for illustration. The approximation v{E,lj) was also used by this 
calculations, namely, the equations (3). It is seen that in the discussed case 
when ± 0, T n = T n o ^ const , the temperature T e (E,u > ) is reaching its 
stationary value in about 50 seconds. However, by solution of the full sys- 
tems of these equations, when = 0, the temperature is growing up all the 
time and almost linearly when t > 1000 sec. 

A system of four equations (22), (23) characterize a full ionized magne- 
toplasma, namely: N e — N , , u ei ^ 0, N n ~ 0, u en ~ v in ~ 0. This is the 
case closely realized in the magnetosphere. However, this system of equa- 
tions has not a stationary solution too. It happened because the loss of the 
energy of the electrons accelerated by the electric field is equal to the energy 
transmitted by them to the ions. The determinant of this system A 2 = 0. 
The temperatures T e and T, are growing up all the time. 

Thus, the limits of applicability of the microscopic approximation used in 
many studies and here, by considering the behavior of all the temperatures 
T e , Ti, and T n by the two approximations: v = const and u(E,uj) should be 
known. Especially it was important to search in this study the role of the 
neutral particles. The results of appropriate calculations are illustrated in 
the next section. 

The heating of the magnetoplasma is accompanied and regulated, in ad- 
dition to the collisions i^i, v en and o in , by many other processes. Their role 
should be of different degree. By a complete solution of this problem, they 
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should be taken into account. The general solution of this task will be- 
come stationary and the growth of the temperature can be stopped. How- 
ever, it becomes a very complicated theoretical task. But, it seems that 
one of these most important processes is the ionization of the neutral parti- 
cles by the accelerated electrons under the action of the electric field E 0 e lwt . 
Other losses of energy which can stop the growth of the temperature of the 
magnetoplasma are, for example: 

the thermal emission of the constituent particles of the plasma i.e. the 
volume Stefan- Boltzman emission of the plasma, 

the heating kinetic losses of the plasma, i.e. the electron-electron and 
ion - ion collisions v tt and z/ u and possibly other kind of kinetic thermal 
losses which, in a sense, are similar to the Landau damping of e.m. waves in 
a plasma. Let us here only look at the ionization process of plasma by the 
accelerated electrons. 


The velocity of the electrons becomes in the ionosphere, depending on 
the frequency u> and altitude Z, larger than the ionization potential E{ of 
the neutral particles. Namely, the thermal, chaotic velocity u e (E) plays the 
prime role by this process. It is considerably larger than the directed velocity 
V e (E). For instance, in an isotropic plasma 


v e (E) = 


2«Te(-£’o, a;) 


1/2 
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iMK(E)} 


20R e {V e (E) 


(38) 


The ionization potential of the atomic hydrogen Hi particles (they are the 
main constituent at the high altitudes of the ionosphere) E, = 13.54 eV. 
Therefore, the process of the ionization begins when v e (E) > 2.2 • 10 8 cm/s. 
The additional density of the electrons A ,r e (E) can become even considerably 
larger than the regular density of the plasma N e when the amplitude of the 
electric field is rather small, even when E 0 = (1 — 2) mV/m. 
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Indeed, the cross-section <7 e (E,) of ionization of the atomic hydrogen in- 
creases rapidly with E, up to a sharp maximum <7 e>mai ~ (7 — 8)10 -17 cm 2 
and then it is smoothly decreasing up to cr etS(a( ~ 2 • 10 -17 cm 2 (Massy at al. 
1969 [11])- At Z = 300 km, for example, where the atomic hydrogen density 
N n ~ 3 • 10 9 , the production of the electrons 

I t = ff e • N n • t> e (E) > (10 to 50) s _1 (39) 


The additional ionization is described by the recombination equation 


dN e 

dt 


LN f - a £ Nl 


(40) 


Growing up in time N e (t) reaches a maximum by = 0 equal to N = 10 9 , 
even if the coefficient of recombination ot e ~ 10 -8 cm 3 • s -1 . However, at the 
discussed here altitudes Z, a e < , <C 10 -8 . Thus, under the action of the elec- 
tric field, the surrounding plasma can become full ionized. As a result, the 
collision frequency between the electrons and ions ME) oc A^ e (E)/T e (E) 3/2 
becomes very large. The growth of T e is stopped very quickly because the 
electron concentration is very quickly increasing. Indeed, in the point of in- 
flection of the time dependence of the iV e (E, f), where = 0, N e = 
Thus, at that point is very large (the value of a e is very small) 

and the tangential of N e (E,t) is almost vertical. 


The discussion of all the noted above points in detail is beyond the frame 
of this study. It is even impossible to do this because many initial data are 
not known and may be found only by treatment of appropriate experimental 
data. It is also seen below that, in many cases, they are not of the decisive 
importance in the altitude region of the Active mission. 
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III. 2 Altitude and frequency dependencies of the temperature 

in the ionosphere. 


The numerical calculations of the temperatures were done for an iono- 
sphere model given in Section II (see page 6). The frequency band was 
widened, namely, F = (1 to 10 4 ) Hz and more because of the general 
interest to study the heating of the ionosphere and magnetosphere by ELF 
waves 0 < F < Fb , Fg is the ion gyrofrequency. By the preliminary calcu- 
lations, it was learned that the growth of the temperatures T e and T, with 
altitude is very large at frequencies F = (1 to 100) Hz and becomes about 
(10 2 — 10 3 ) and many times larger than the temperature of the initial plasma 
T n0 . Therefore, F = 1 Hz is used to characterize this band of frequencies. 
It illustrates more dramatically the region where it is necessary to stop in 
the theory the temperature growth by including additional losses of energy. 
The zone of ionization by the accelerated electrons is marked on the plots 
by a line and vertical arrows. Besides, it was learned that in frequency band 
F > (10 2 tolO 4 ) Hz and more, the altitude dependencies of the temperature 
are qualitatively the same, therefore, the frequency F = 1000 Hz was used 
to characterize this frequency band. The altitude temperature dependence 
on this frequency is sufficiently typical for this frequency band. These two 
frequencies are used also to characterize the angle 0 dependence of the tem- 
peratures. Let us remind here that © is the angle between the electric E and 
magnetic Bq fields. 


a. Numerical results. Approximations v = const, v — v{E,u>). 

Altitude dependencies T e /T n0 and T,/T n0 of the electron and ion tempera- 
tures for 0 = 0, F = 1 Hz and F = 1000 Hz are shown on Fig. 12. They were 
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calculated by the two approximations v = const and v e = t'(T) = i/(E,w) 
(see (24)) by the selfconsistent solution of the system of equations (22) and 
shortened system of equations (23) when ^jf- = 0, T n = T no . Altitude 
dependencies by approximation v = const, ^2*. = 0, T n = T n o are also 
shown: for 0 = 0 and four different frequencies F = (1 and 100) Hz, 
F = (10 3 and 10 4 ) Hz on Fig. 13, for F = 1 Hz and angles 0 = 
(0, 30, 66, 75, 85 90)° on Fig.14, and for / = 10 3 Hz and angles 0 = 0°, 
0 = 75° calculated by the approximations v = const and v = v{T) = v(E,u>) 
on Fig.15. 

By examination of all these figures, it is seen the following: 

1. The character of the altitude dependencies of the temperatures is differ- 
ent in the frequency bands F = (1 to 100) Hz and F = (10 3 to 10 4 ) Hz. 
This is also seen from Fig. 16, where frequency dependencies are shown 
at Z = 300 km for the approximations v = const , v = i/(E,u>) when 
^ = 0, T n = T n o. 

2. The temperatures of the electrons T e and ions T, are growing up quickly 
on both frequencies. They become equal in the altitude region Z ~ 
(200 — 250) km. In this region ( T e and 7)) > (10 2 — 10 3 )T n o. 

3. The quick growing of temperatures, when F = 10 3 Hz is slowed down 
and at Z ~ 200 km they have maxima by the approximation v = 
v{Eq,uj) and v — const. 

4. At Z > 300 km the temperatures are smaller by the approximation 
v = v{Eq,u> ) than by the approximation v — const. 

5. The process of ionization by the accelerated electrons is acting almost 
in all the altitude region when F — 1 Hz and does not play a role, 
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Fig 12. Altitude dependencies of T e /T n0 and T{/T n 0 , 0 = 0. 

Approximation v = const, F = (1 and 10 3 ) Hz — * 1 0 , 2 0 . 
Approximation v = y(E,u;), F = (1 and 10 3 ) Hz — > 1e, 2 e- 
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Fig 13. The same as on Fig. 12 T e /T n0 . Approximation v = const. 
Frequencies F = (0, 1 0 2 , 1 0 3 , lO 4 ) Hz — + 1., 2., 3., 4. 
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Fig 14 . The same as on Fig.13. Angles between E and Bo: 
0 = (0, 30, 60, 75, 85, 90)° -» 1, 2, 3, 4, 5, 6. 
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Fig 15 . The same as on Fig. 12, F = 10 3 Hz. 

Approximation v = const, 1 0 — » (0 = 0), 2o - ♦ © = 75°. 
Approximation v = i/(E,u>), ljr — * (© = 0 )' 2e — + 0 = 75 . 
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10 100 1000 100C 


F , Hz 

Fig 16 . Frequency dependencies of T e (E,u>)/T n0 , Z = 300 km 
Approximation v = const, © = (0, 75, 85)° — * 1., 2., 3. 
Approximation v = ^(E,u)), 0 = (0, 75, 85)° — > l#, 2 e, 3,e. 
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when F > 10 3 Hz, especially at Z > 300 km. This process, namely 
T e ( E,u») > l0 2 T nO could be stopped by the theoretical calculations 
when F = 1 Hz already at the altitude Z = (120 — 150) km. This 
does not mean that the temperatures of the electrons and ions can not 
become larger than the critical temperature T e ,(») of ionization, let us 
say, of the atomic hydrogen#!, which is equal to 

T et(i ) = 13.5 eV = 1.16 • 10 4 • 13.5 = 1.57 • 10 5 , K°. (41) 

To answer this question, theoretical studies must be done, taking into 
account losses of energy, in particular, noted in Section III.l. 

6. The influence of the magnetic field becomes effective by angles 0 > 30°. 
Thus the given above formulae for an isotropic plasma may be widely 
used for evaluation of the temperatures T e and T; in a magnetoplasma. 


b. ^=0and^*0. 

The calculations given in Section Ilia are in general confirmed by ex- 
amination of the results of calculations of the temperature by selfconsistent 
solution of the full systems of equations (23) and (24), i.e. when ± 0, 
T n = T„(E,uj) 4- T n0 . These results of calculations are illustrated by Fig. 17. 
Not to overload this figure, only the T e and T n dependencies are given in it. 
The main new important characteristics of these dependencies, in addition 
to the given above 6 points, are the following: 

1. The temperature T„ of the neutral particles is growing up very quickly 
at altitudes Z > 200 km. It approaches the temperature T e of the 
electrons at Z > (400 — 500) km and becomes close to T e similarly to 
the temperature of ions, 7 1 ,. 
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Fig 17. Altitude dependencies of the temperatures. 
Approximations v = 0 = 0. 

T e (E,u)/T n0 : F = 1 Hz, = lr - ^ # 0, 

F=l0 3 Hz, 2. -»^ = 0, 2 r -^^0, 

F = 1 Hz, 3. -> = 0, 3 r -> ^ # 0. 
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2. At the altitudes Z > 400 km, the temperatures T e are larger on both 
frequencies by the approximation ^ / 0 than by the approximation 

= 0 . 

3. At least at Z > 300 km, the heating of the neutral particles becomes 
a source of the heating of the electrons and ions. The growth of T n 
should be stopped, similarly to the case = 0 at these altitudes. 

4. In many cases, the heating of the neutral part, ’es is an important part 
in the process of heating of the magnetoplasma. 


Summary 

Results of calculations of the moduli of the electric field \E\ generated by 
an electric dipole at different altitudes of the ionosphere Z ~ (300 to 2500) km 
up to the bottom of the magnetosphere Z ~ 6000 km and in the frequency 
range Fb <C F < fs are given. It is supposed that the source of this field 
is placed on a satellite, moving on an elliptical orbit which crosses this alti- 
tude region of the magnetoplasma. The electric field is recorded on another 
Sub-satellite moving around the source at distances of (10-s to ~ 100 Arm) 
or a little more. The linear theory for a homogeneous medium may be used 
for treatment of the experimental data. The altitude frequency and an- 
gle dependencies of the field are studied in all the regions of their largest 
enhancement. They are the Axis field |£o|) jd ~ 0, close to the direction 
of the Earth’s magnetic field line Bo, the fields |£sf| and lEfiesStl i n the 
Storey cones, ft ~ (0 to 10 — 20) degrees, the Resonance cone field \Erc 3 \, 
~ (0 to 90) degrees ; /? is the ray direction angle. The Axis field |Eo| is 
studied here in detail. It is the most interesting and important peculiarity of 
this problem. The fields |£st|) I EjievSt 1 , and j Eft es | are increased by a factor 
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of (10 to 10 5 ) and more by approaching the low hybrid frequency. The value 
of |2?o| is increasing at the altitude (500 — 2500) km about 10 6 times and it is 
concentrated in angles A/? < 10 -1 degrees. These waves produce formations 
similar to laser beams. In the used model of the ionosphere, it should happen 
in our case in the frequency range F ~ (20 to 10) kHz. 

The extended recently microscopic theory of heating of the magneto- 
plasma is used to calculate the temperatures in the ionosphere. Adequate 
formulae and some numerical results are given, particularly, based on the self- 
consistent solution of a system of equations of the derivatives velocities dV J t | l | n 
and of the temperatures — All kind of collision frequencies between the 
constituents of the magnetoplasma v et , u en and are taken into account. 
Numerical results illustrate the altitude and frequency dependencies of the 
temperature in the frequency band F = (1 — 10 4 ) Hz and altitude range 
(100 to 10 3 ) km. The heating of the ionosphere is growing up very quickly, 
with altitude. The temperatures reach the ionization potential even when 
the electric field \E\ is about some units of mv/m. At Z ~ (150 — 200) km, 
T e > 10 2 T nO in this region. 

It is very important to investigate in detail the studied in this report 
problems by experiments on satellites. The treatment of the experimental 
data by the given above results of calculations will help both, to understand 
more the physics of different processes in the magnetosphere and also to 
improve the existing theories of these phenomena. 
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